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We present a systematic way to combine the global color model and the 
instanton liquid model to calculate the chiral Lagrangian coefficients. Our 
numerical results are in agreement well with the experimental values. 
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Quantum chromodynamics (QCD) is successful to describe strong inter- 
action in terms of basic force at high energy. However, the large coupling 
constant at low energy makes the perturbation theory break down. At low- 
energy scale, the freedoms of quarks and gluons realize themselves in hadronic 
level as mesons and baryons, which is characterized by chiral perturbation 
theory (>PT). [1 ' 2] The xPT can give a lot of satisfactory explanations for 
experimental data while the coefficients of the chiral expansion are input 
parameters and their values are determined from experimental observables. 
Since the present technique cannot derivate the chiral Lagrangian from QCD 
directly, a class of chiral invariant QCD based models are used in deriving 
the chiral Lagrangian which is characterized by the effective quark-quark 
interaction. 

The global colour model (GCM), which has provided a lot of successful 
descriptions of the long distance properties of strong interaction and the QCD 
vacuum as well as hadronic phenomena at low energy, is based on a trun- 
cation of QCD Lagragian preserves global colour gauge transformation and 
permits a large N c expansion.^' ^' ' ^ The instanton liquid model (ILM) 
of the QCD vacuum is based on a semiclassical approximation, in which all 
gauge configurations are replaced by an ensemble of topological nontrivial 
fields, instantons and anti-instantons.''' 7 ' ^ Although it does not give rise to 
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a long range confining force between quarks, the instanton vacuum has been 
shown to provide a good phenomenological description of many hadronic 
properties. Therefore it is interesting to calculate the chiral coefficients with 
the two coupling nonperturbative models. In fact, the method of the two cou- 
pling nonperturbative models to calculate physical quantities already exists 
in the literature and is successful in describing some hadronic phenomena. 
In this Letter, we approximate the effective quark-quark interaction by a di- 
lute instanton liquid ensemble and calculate the chiral coefficients from the 
GCM following Ref. [4], while the pioneer work for obtaining chiral hadronic 
phenomena through GCM was given by Roberts et a/J^ 

Here we describe briefly the main skeleton of the GCM. The global colour 
model, based upon an effective quark-quark interaction which is approxi- 
mated by the dilute instanton liquid configurations, is defined by a truncation 
of QCD Lagrangian while maintaining all global symmetries. 

The generation function for QCD in Euclidean space can be written as 

Z[fj,rj] = J VqDqVA exp j— S + J d 4 x(fjq + , (1) 



where 



S = Jd 4 x{q ^^-ig^A^+if; 




and G a = d^A^ — d u A a ^ + gf abc A b ^Al. We leave the gauge fixing term, the 
ghost field term and its integration measure to be understood. Here A^ are 
approximated by the fields of the dilute instanton liquid. In order to preserve 
the chiral invariance of the above action S (not the full action), the quarks 
have to be coupled to a external field ip(x) in stead of current quark mass m, 
the field ip(x) transforms in a certain way under chiral transformations^' ^> 4] 
(For detailed discussion one can see Ref. [4]). 
Here we introduce 

e W[J] = J VAe Jd±x(-l G ^G* v +J*A*)^ ( 2 ) 

with J^(x) = q(x) , y^q(x) and rewrite the generation function as 

Zty, 77, 77] = f VqVqe~ f d4x W Td+ ^ q - m -^ e w[J}^ (3) 
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then expand W[J] in the current J^: 

W[J\ = \J d 4 xd 4 yJ«(x)D« b u (x, y)J b M + ± / J^D^ + • • • . (4) 

The GCM is defined by a truncation of the functional W[J] in which the 
higher order n(> 3)-point functions are neglected, and only the gluon 2-point 
function D^ u (x, y) is retained. This model maintains all global symmetries of 
QCD and permits a 1/N C expansion, however local SU(3) gauge invariance is 
lost by the truncation. For simplicity, we use a Feynman-like gauge D ab v {x — 
y) = b[ivb ab D(x — y). Using the functional integral approach, we introduce 
bilocal field integration variables B d (x,y) and integrate out the quarks and 
gluons fields to change dynamical freedoms of fundamental QCD to hadronic 
variables. This method has an obvious advantage that the effective hadron 
field interactions retain knowledge of their subhadronic origin. The resulting 
function Z GCM takes the following form: 

Z GC u = J VBe~ s W, (5) 

where the action is given by 

S[fj,r,,B\ = -TrLniG-^ J ^y ^^f^K J d*xd 4 yfj(x)G(x,y; [B e ])r](y), 



(6) 



2gW(x-y) 

and G' 1 is defined as 

G-\x, y) = ( a .d + ^{x))5{x -y) + A 9 B 9 (x, y). (7) 



Here the quantity A arises from the Fierz reordering of the current-current 
interaction term in Eq. (4) and represents the direct product of Dirac, flavor 
SU(3) and color matrices. 

The vacuum configurations are defined by minimizing the bilocal action: 
= 0, following this mean field approximation, we obtain 

?7,T)=0 

B e (x -y)= g 2 D(x - y)tr[A 9 G (x - y)\. (8) 

These configurations provide self-energy dressing of the quarks through the 
definition E(p) = A e B e (p) = 27 • p[A(p 2 ) - 1] + B(p 2 ). In fact, Eq. (8) is the 
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rain-bow Schwinger-Dyson equation. In this Letter, the Schwinger-Dyson 
functions A(p 2 ) and B(p 2 ) are approximated by the results of ILM through 
equating the two point Green functions of the corresponding quark in the 
two nonperturbative methods. 

With the GCM generation function Eq.(5), it is straightforward to cal- 
culate the vacuum expectation value of any quark operator of the form 
O n = (gi 1 A^J- 1 q 7 - 1 )(g i2 AQ 2 gj 2 ) • • • (ft„AQ n ^J in the mean field vacuum. Here 
AW represents the operator in Dirac, flavor and color space. 

Taking the appropriate number of derivatives with respect to the external 
sources rji and fjj of Eq. (5) and putting rji — fjj — , we obtain 

<0| : O n : |0> = (-)"E[(-) P AS 1 • • • Atl(Go) jlipW • • • (Gok M J, (9) 
v 

where p represents the permutation of the n indices. In fact, the above 
equation can only be correct with the assumption of vacuum saturate or 
vacuum dominance. ' ^ Furthermore, in calculation, we should subtract the 
perturbative part of the two point Green function for the quark,' 1 which is 
here neglected due to the small current quark mass. 

The quark condensate (0| : q(0)q(0) : |0) can be inferred from Eq.(9), 

(0|:^0MQ ): [0) = -^f^ ^W B2(s) . (10) 

In the instanton liquid model, A(q) = 1. Here /x is taken to be 1 GeV 2 , 
the scale dependence of the quark condensate is thoroughly discussed in the 
previous paper.' ' The quark condensate is calculated for later use in the 
expressions for the Chiral coefficients. 

We expand the bilocal fields about the saddle point, 

B 9 (x,y)=B e (x-y)+B e (x,y), (11) 

and then examine the effective interactions of the fluctuations, B. The bilo- 
cal fields contain information about internal excitions of the qq in addition to 
center of mass motion which is associated with the usual local field variables. 
Thus, they can be projected onto a complete set of internal excitations Y e n 
with the remaining center of mass degree freedom represented by the coeffi- 
cients 4> e n {P). In fact, the coefficients 4> e n {P) can be identified with the meson 
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fields. The bilocal fluctuations can be written as 

B e (P,q)=T,^n(PW 9 n(P,<l)- (12) 

n 

The functions T 6 n are in general eigenfunctions of the free kinetic operator 
of the bilocal fields. At the mass shell point, P 2 = —M 2 , they satisfy the 
homogeneous Bethe-Salpeter equation in the ladder approximation for the 
given quantum number 9 and modes n.' ' 

To obtain the chiral coefficients, we take a digression to discuss Goldstone 
theorem. Dynamical chiral symmetry breaking is associated with the occur- 
rence of a massless Goldstone mode which can be related to the dynamical 
generated scalar self-energy of the quark propagator in the limit of vanish- 
ing current quark mass. By analysing the corresponding zero momentum 
behaviour of the axial- vector Ward identity in the chiral limit, we can ob- 
tain the Goldberger-Treiman relation for the quark-pseudoscalar vertex and 
then reduce the bound state function T e n (P, q) from ladder Bethe-Salpeter 
equation for the pseudo-scalar Goldstone mode with P 2 = to the quark 
self-energy function.^ 

Our main goal here is to obtain the effective action for the Goldstone 
mode, the higher mass fluctuations present in the bilocal fields are neglected. 
In fact, the higher mass fluctuations should have some effects if they have 
been integrated out. Based on the above discussion, we can write the full 
bilocal field of quark fields in the following formal : 

A e B e = E(x -y) + B(x - y)[U B fi±V) - 1], (13) 

where U^(x) = PrII (x)+PlU'(x) with Pr^ the standard right-left projection 
operators. For the SU (3) flavor case under consideration here the chiral field 
U is defined as U = e iXa W f " . 

Next, we follow Ref. [4] to calculate the chiral coefficients and we con- 
sider only the real contribution to the effective action. The pioneer work for 
obtaining the chiral hadronic phenomenon through global color model was 
given by Roberts et al.^' The restriction of the fluctuations to the Goldstone 
modes with UW = 1 entails that the second term of the action in Eq.(6) is 
independent of the filed U and can be neglected. Thus, the real contribution 
to the action can be written as 

S = Re[S] = -l-TrLniG-^G- 1 }^), (14) 
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where 



G- 1 = 7 • d x A(x -y) + iK^)6(x -y) + B(x - y)U b {^-). (15) 



Here A(x — y) — 5(x — y) in the ILM. 

By expanding the logarithm and dropping the irrelevant constant terms, 
we obtain the chiral Lagrangian to the fourth order (in Euclidean space): 



where x( x ) = ~~ ^(qq)^^)/ fn an d the remaining trace is over flavor. The 
coefficients L\, L2 {Li = 2Li), L3, L5, L§ and it decay constant f n are func- 
tions of the scalar part of quark self energy and its derives, as the expressions 
are lengthy, one can see Ref. [4] for details. 

In performing the gradient expansion it is important to keep the x- 
dependence of the external field ip{x). After carrying out the gradient expan- 
sion to the fourth order, we employ the equation of motion which is obtained 
at the second order and depends on the external field ip(x). Finally we can 
identify ip{x) with the current quark mass matrix. Failure to keep the x- 
dependence of ip to the very end violates chiral invariance and will render 
some unphysical results for the low-energy coefficients. 

From the above expressions, we can see that the coefficients L4, Lq, Li 
vanish, which is consistent with the large N c QCD result. Some modifications 
may be obtained if we have integrated out higher mass states. 

Instantons, as the solutions of the classical Yang-Mills equation of motion, 
give many successful descriptions of strong interaction phenomena, such as 
U(1)a problem, dynamical chiral symmetry breaking, tunneling the vac- 
uum and so on. After the great achievement of 't Hooft's computing the func- 
tional determinant in the one-instanton background, a tempt was made to 
understand quantitatively the nonperturbative QCD vacuum configurations 
in terms of instanton ensembles.' ' ' However, the notorious infrared prob- 
lem due to large size instanton made many researchers no longer believe in the 
dilute instanton gas model. To circumvent the disaster, one may think the 



/ d 4 x{ftr[d»Ud^} - tf[U X * + XU*} 

-L^tr^Ud^]) 2 - L 2 tr[d^Ud p U^tr[d^Ud v U^ 
-L z tr[^Ud^d v Ud v U^\ + L 5 tr[d^Ud^{U X ^ + X^)} 
-L 8 tr[xU^xU^ + Ux^Ux^}}, 




s = 



(16) 
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overcome that the vacuum gets more filled with instantons of increasing size 
for larger distances, at some scale there might be some repulsive interactions 
to stabilize the medium while the semiclassical treatment is still possible and 
the instantons are not much deformed through their interactions, thus form a 
dilute instanton liquid ensemble. In Ref. [13, 14], the authors used variational 
techniques and mean field approximation to deal with the statistical mechan- 
ics of the instanton liquid while stabilizing it by phenomenological hard core 
and repulsive interactions derived from ansatz of gauge field. Subsequently, 
the propagation of light quarks in the dilute instanton liquid ensemble was 
formulated quantitatively and refined due to the treatment of the nonzero 
mode propagation.^' ' Although it does not give rise to a long range con- 
fining force between quarks, the instanton vacuum has been shown to provide 
a good phenomenological description of many hadronic properties.^' ^ Fur- 
thermore, lattice calculations support that the vacuum gluon configurations 
can be approximated by the dilute instanton liquid. ' ^ Its shortcoming will 
be that the vector self-energy function of quark propagator in the instanton 
ensembles A(q 2 ) = 1 which is not consistent with the corresponding one in 
lattice QCD and other potential models. 8] j n addition, the q 2 > lGeV 2 
behaviour of the scalar self energy function B(q 2 ) of quark propagator in 
the instanton medium falls as 1/p 6 which disagrees with QCD. Thus, it is 
interesting to substitute the quark propagator which is dressed through prop- 
agating in the dilute instanton medium into the expressions for the n decay 
constant and chiral coefficients. Our numerical results fit experimental values 
well. The quark propagator in the ILM can be written as 

G(p) = [i 7 • p + Mip)}- 1 + 0(N/VN c ), (17) 

where 

N P 2 <p 2 (p) 

<p(k) = irp 2 ^-(I (z)K (z) - / 1 (^)K 1 (^))| 2=fcp/2 , (18) 

where / and K are the modified Bessel functions, N/V (200MeV)^ 4: and 
p « {mQMeV)- 1 . 

In conlusion, we have provided a systematic way to combine the GCM 
with the ILM.' 1 Then based on the nonperturbative quark propagator de- 
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rived in the ILM, we calculate chiral coefficients which are expressed in terms 
of the scalar part of quark self energy and its derives. 

Our results are compatible with the experimental data. Our results are 
U = S6MeV, (qq) = -(217MeV) 3 , L x = 0.8 x 10" 3 , L 2 = 1.6 x 1(T 3 , 
L 3 = -3.9 x 10" 3 , L 5 = 2.0 x 10~ 3 , L 8 = 1.0 x 10~ 3 . The experimental values 
are L x = 0.9 ± 0.3 x 10~ 3 , L 2 = 1.7 ± 0.7 x 10~ 3 , L 3 = -4.4 ± 2.5 x 10" 3 , 
L 5 = 2.2 ± 0.5 x 10~ 3 , L 8 = 1.1 ± 0.3 x 10~ 3 

Following the GCM, we truncate the QCD Lagragian to include only two 
point gluon Green function which is approximated by the dilute instanton 
liquid ensemble and describe the mesons as fluctuations about the saddle 
point of the effective action. The chiral coefficients obtained in this way 
are consistent with the large N c QCD, and embody their underlying theory 
through the dependence on the quark self-energy and its derives. Some mod- 
ification will be made if we have integrated out the high mass fluctuations 
instead of neglecting them. 

To obtain the nonperturbative quark propagator, we have to resolve the 
rain-bow Schwinger-Dyson equation. However, the form of the gluon propa- 
gator g 2 G{p) in the infrared region is unknown, one often uses model forms as 
input in the previous studies of the rainbow Schwinger-Dyson equation.^' ^> 5] 
Obviously, this dressing comprises the notation of constituent quarks by pro- 
viding a mass m(p 2 ) = B(p 2 )/A(p 2 ) with (m(0) = 418MeV) in the ILM, 
which is corresponding to dynamical symmetry breaking. Here we use the 
dilute instanton liquid to approximate the QCD vacuum configurations, the 
quark propagator is dressed through propagating in the dilute instanton liq- 
uid ensemble and contains valuable nonperturbative information which man- 
ifest itself in the coefficients. Although ILM has shortcoming as mentioned 
above, it does give satisfactory chiral coefficients and can be extended to 
other systems. Furthermore, the lattice calculations support that the dilute 
instanton liquid ensemble is a good approximation of vacuum configurations. 
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